NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



QUARTERLY REPORT NO. 3 


TRANSIENT AND DIFFUSION ANALYSIS OF Hg Cd Te 


1.0 INTRODUCTION 


This is the third report of technical progress on Contract No. 
NAS8-33698, titled, "Transient and Diffusion Analysis of Hg Cd Te". This 
report concerns technical work performed between 20 March and 20 June 1980, 
and includes work being carried out under the increased Scope of Work, 
"Exhibit B" of the modified contract. 

The work during this period has concentrated on the solution of 
the one dimensional planar interface solidification problem. This problem 
has been formulated in a most general way and numerical solutions have 
been obtained. In addition, progress on the design and construction of 
a high gradient furnace will be reported. 

It is useful at this point to review why the one dimensional solid- 
ification problem is of interest. Previous work by Tiller and Jackson (1) 
and Smith, Tiller, and Rutter (2) addressed the problem of solute distri- 
bution in a dilute alloy system with constant segregation coefficient 
and a constant growth rate. By treating the one dimensional planar inter- 
face case in three parts, they were able to obtain analytical solutions 
for solute distribution in the initial transient, steady state region, 
and final transient. As pointed out in the previous progress report. 


these solutions do not apply to the Hg Cd Te system due to the following 
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behavior that violates most of the previous works assumptions. In 
contrast to the dilute alloy system, the Hg Cd Te system has for a 
initial liquid solute mole fraction of . 2(Cd Te) : 

• A solute mole fraction in the solid that varies 
from about .6 for the first to freeze to .2 at 
steady state. 

• A changing interface solidification temperature 
from 8Q0°C Initially to 705°C at steady state. 

• Widely separated liquidus and solidus curves with 
a segregation coefficient k that varies from 2.9 
to more than 6 over the approach to steady state 
growth . 

The goal here is to reformulate the problem to allow for changing 
k and interface temperature. Recent work by Favier (3) has pointed 
out that the interface velocity cannot remain constant throughout a 
solidification experiment and that the applied thermal field in effect 
determines the solidification rate. Favier has solved the one dimen- 
sional problem for a rather complicated thermal field but assumed a 
constant value for the segregation coefficient, thus limiting the 
applicability of this work to the Hg Cd Te system. In the following 
sections we will formulate the general problem, indicate a numerical 
method for solution of the problem and apply the results to experimental 
solute profiles to determine the effective diffusion constant for 


2.2 TRANSIENT ANALYSIS 


2.1 ANALYSIS OF THE PLANAR INTERFACE CASE 

We wish to formulate the most general one-dimensional solid- 
ification problem in keeping with the Bridgam - Stockbarger growth 
configuration and the HgCdTe solid solution system. We shall make the 
following assumptions: 

I. The thermal and concentration fields considered are one- 
dimensional. 

II. No significant diffusion occurs in the solid. 

I 

III. The thermal field is a specified function of the spatial 
variables and time and does not change as a result of the 
solidification process . The thermal field is taken to be 
monotonic in the spacial variable. 

IV. The sample movement rate, R, is constant. The growth rate, 
v, is not necessarily constant. 

V. We will deal with a binary or pseudobinary system. 

We shall choose spacial variables such that z is measured from a 
fixed laboratory frame and x is measured in the sample from the first 
to freeze point (ampoule tip) . The variable q is the distance into the 
liquid from the solid - liquid interface (see Figure 1.). 
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The starred (*) values refer to variables evaluated at the solid- 

liquid interface. Choose t * 0 to be the start of solidification in 

& 

the tip of the ampoule and z (0) the z - coordinate of the ampoule 
tip at t * 0. We may then relate the various coordinate systems by 


X * z 

+ Rt 

- z* (0) , 

(1) 

* 

X * 

* 

z + 

Rt - z*(0) , 

(2) 

and 


* * 

(3) 

n ■ x 

- Rt 

+ z (0) - z (t) . 


The total length of the sample is SL, and n * x denotes the 

end of the ampoule in the n system. 

Let C(n,t) denote the solute mole fraction in the liquid ahead of 
the interface. It is more convenient to express the thermal field in 
the fixed coordinate system, thus, 

T = T(z, t) (4) 

and we shall assume a thermal field of the form 

T * T + Gz (5) 

o 

where G is: thus the thermal gradient, assumed to be constant over the 
region of interest. 

We define the liquidus and solidus components of the phase diagram 
(binary or pseudobinary) by L(C) and S(C) respectively. We have 


S(0) = L(0) and S(l) = L(l) . 


( 6 ) 



Furthermore, we may define the solid - liquid interface location 

* - V 

z (t) such that 

T(z\t) = L(C(0, t) ) (7) 


since n * 0 denotes the interface location in the n system. 

Writing the diffusion equation in the z - coordinate system and 
transforming into the n system gives : 



with the associated boundary conditions 



(9) 

( 10 ) 

(ID 


The boundary condition in equation (10) is normally written in a 
slightly different form. If we consider a semi-infinite sample, 
equation (10) may be written as 


C(n»t) « C as n -*■ 00 for all t. (12) 

This form of boundary condition is sufficient for the semi-infinite 
sample and allows computation of the initial solidification transient 
and approach to steady state growth. It is our aim here to solve for 
both the initial and final transients for a finite length sample, hence 
requiring that the first derivitive of the concentration in the l'iquid 
vanish at the far end of the ampoule is necessary to obtain the initial 
transient, steady state region, and the final transient. Equations 


(10) and (12) are obviously equivalent as JL->- • . 

The boundary condition in equation (11) arises from conservation 
of total solute and is applied at the interface n“ 0. Here k(C) is the 
segregation coefficient 


C (0,t) 

k(c(0 ’ t)) * C (0, t) 

la 


(13) 


where S and L denote the solid and liquid sides of the interface re- 
spectively. Note that we have not assumed a constant segregation co- 
efficient, but k depends on C(0,t) and is defined by the phase diagram 
of the system under study. The growth velocity v(t) is related to the 
sample pull rate R by 

* 

v(t) « R + , (14) 

* 

where we recall that z is defined by equation (7) . 

In addition to the equations above we shall assume that the system 
under study is such that we may approximate the segregation coefficient, 
liquidus portion of the phase diagram, and thermal field by the _ 
following: 


k * 

a + £5C + 5C 2 

(15) 

L = 

T a + m a c + yc 2 
A A 

(16) 

T * 

T + Gz 

a 

(17) 


Equations (16) and (17) now allow one to solve for via equation 

(7) 

T(z*,t) = L(C(0,t)), 
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giving 


z*(t) - (1/G)(T - T + M.C(0,t) + Y C 2 (0,t)) 

A O A 


and 


- (1/G)(M a + 2 YC(0,t)) 5 - C( °^ - 


( 18 ) 


In summary, under conditions outlined above we seek to solve 

D 32c(n - ; - t) + (R + |f- ) Maill , . n . > o (19) 

. z at on ot 

on 

subject to 


C(n»0) = C q , (20) 

jCfit fc > . = 0 at n * &* for all t > 0 , (21) 

on 

(R + ||- )(a + 0C(O,t) + 6 C 2 (0,t) - l)C(0,t) - D - C( ^ - 
at n " 0 . (22) 


The numerical approach to the solution is oulined in the next section. 
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2.2 NUMERICAL APPROACH 


The differential equation and boundary conditions outlined 
in the previous section can be viewed as an initial value problem 

t 

in time and a boundary value problem in the spatial variable. This 
set of equations may be solved by finite difference techniques, but 
we consider a different approach to be more appropriate. The al- 
gorithm chosen will be outlined below. 

Let the incremental time and spatial steps be At and An tespect- 
ively. We will assume that at some time t^, the values of C(n.t) 

have been determined for tit and all n. The goal is to then com- 

o 

pute C(n.t') , where t' * t + At for all n. Assume a value for the 
concentration at the interface at t * t’, say 

C(0,t') = C* , (23) 

* 

where C is a trial value. We then compute 


9C(0,t’) C - C(0,t ) 

o 


3t At 

This allows us to determine the interface velocity v by 
* 


ff- - (l/G)(M A + 2 ye (0 , t ' ) ) 


and 


v * R + 


dz 

dt 


(24) 


(25) 


(26) 


The first derivitive in n may be found from the boundary condition 
(Equation (22)) as 

A 

-■-■ fe -- - “ (1/D) (R + —-)(« + BC<0,f) + 6C 2 (0,t , )-l)C(0,t’)) 


The second derivitive of C in n is found via the diffusion equation 
(Equation (19)) giving 

aic . tQjt 1 ) « (1/D) ( S&gJL'J . . ( R + ||_) ) ( 28 ) 

3n 

We then may use a Taylor series expansion to obtain our desired results. 


c(An,t t ) - c(0,t') + + l- c .CQi. t,) An_ 

8n 3n 2 2 


(29) 


as the value for the concentration at q*An in front of the interface 
at t - t'. 

So far we have simply computed C(n,t) at the first spatial in- 
crement away from the interface. We now wish to propagate this value 
out in n until the end of the ampoule is reached (n * £') and test 
to see if the boundary condition 


dC , | N 

■ 3? <"•* > 


0 at ri= 


(30) 


is met. The solution may be propagated to n' * 2An, the next spatial 
step, by first calculating 


3C(An,t') ^ c(An,t') - c(An,t Q ) 


(31) 


at At 

We then may use a Taylor series expansion and the diffusion equation to 
obtain 


3C 

3n 


(Ay. t*) . ^ t') + g CC0,t'. ) 
Un,C } 3n ' 2 . 


An 


(32) 


an 


(An ft') = (l/D) ( || (An.f) - (r +|5-) — ) . (33) 


We note here that the value of the Interface velocity R + -r-— does 

at 

not change as we increment the spatial variable. A second expan- 
sion leads to the desired results 

C(n',t’) - e(An,t') + An+ • 9 - 2 ^< Ar lr t, . ) . ) (34) 

3n 

The process from equation (31) to equation (34) is repeated for the 

next spatial step 3Aq, act. until ti* £'» the end of the ampoule at 

this time step. A test is made to see if the boundary condition 

* 

(Equation (30)) is satisfied. If not, a new trial value of C is 

selected and the entire process involving equations (23) to equation 

(34) is repeated until the boundary condition is satisfied. 

We began this explanation by assuming all of the C(r),t) were 

known at some t ■ t . The initial value boundary condition, 

o 

C(rj.O) = C , allows one to take t =0 and proceed in a straight- 
o o 

forward manner to calculate C(n,t) for any t. At each time step, 

* 

after the value of C(0;.t) - C has been found that satisfies the bound- 
ary condition in Equation (30) , we may then proceed to obtain the 
solute concentration in the solid. 

The interface position, measured from the tip of the ampoule, is 
given by x (Equation (21)) where we note that z (t),and thus x, 
is determined by the thermal field and the liquidus curve of the phase 

diagram evaluated at C(0,t). The solute concentration in the solid 
* 

C g (x ) , is defined by 

C (x*) - k(C(0,t))C(0,t) 

S 


( 35 ) 


where k Is assumed Co be known via Equation ( 15) . After the solid 

* 

concentration is determined at a position x , the remaining amount of 
liquid to be solidified is found as 

« 2,- x . (36) 

The time is incremented by At, and the entire processes repeated tc 

* * 

give new values of x (t) s C g (x ), ect., until the end of the ampoule 
2' ■ 0 is reached. In the following section, representative plots 
of these variables and fits to actual experimental data will be presented. 


2.3 RESULTS 


la this section we will present typical results obtained from the 
previous analysis. Figure 2 is a plot of the mole fraction of Cd Te 
in the solid versus the length of the sample measured from the first 

to freeze tip. This curve is calculated assuming a sample pull rate of 

“ -5 2 

.1116 cm/hr for a value of the diffusion constant of 5 x 10 cm /sec. 

The length of the sample was taken to be IS cm. Data necessary to fit 
the liquidus portion of the phase diagram and the segregation coefficient 
was supplied by Dr. S.L. Lehoczky of MacDonald Douglas Research Labora- 
tories in St. Louis, Missouri and his co-workers. The thermal gradient 
was assumed to be a constant 200°C/cm over the region of interest. 

In figure 2 we note the initial transient characterized by decreasing 
Cd Te content as the diffusion layer in the liquid was being built up; 
the steady state region in which the solid is of uniform composition, 
and the final transient in which the solute is depleted in the remaining 
liquid resulting in a decreasing solid solute composition. 

We may now compare calculated curves such as figure 2 to the ex- 
perimentallyddetermined composition profiles obtained by density measure- 
ments or X-ray microprobe analysis. The fit to the liquidus curve and 
segregation coefficient are determined by the plase diagram, so the 
diffusion constant is the only free parameter in analysing the composition 
profile. Figure 3 is the same theoretical data as in figure 2, plotted 
on an expanded split scale. In this case we have included experimental 
composition profile data by Lehoczky, et. al.of MDRL for Hg Cd Te 
crystal grown at a pull rate of .1116 cm/hr. The two theoretical curves 
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-5 2 

are for values of the diffusion constant of 5 x 10 cm /sec and 
-5 2 

7 x 10 cm /sec. Ue note that the agreement between experiment and 

theory is quite good in this case, with the best fit obtained for 
-5 2 

D ■ 5 x 10 cm /sec. The first three experimental points lie well 
off the theoretical curve. The most likely explanation for this is 
that in this region, the crystal is solidifying within the conical 
ampoule tip and thus a one dimensional calculation does not adequately 
reflect initial solidification conditions. We also note that the sensi- 
tivity to the diffusion constant is good enough to allow determination 
of the diffusion constant to approximately 15% if the scatter in the 
experimental data is small. 

One of the primary goals of this study was to analyze the solid- 
ification process in Hg Cd Te so that growth conditions may be optimized. 
Figures 4 and 5 are an example of the value of this analysis in choosing 
proper growth conditions, and are calculated for the same parameters 
as the previous graphs. In the plot of interface velocity versus time, 
we see that it requires almost 25 hours of growth before the growth 

rate is equal to the sample pull rate. Of more practical use is figure 5 

* 

which gives the interface location z (t) as measured in the coordinate 
system fixed to the furnace mid-plane. We see that for the assumed 
mid-plane temperature of 700°G and thermal gradient of 200°C/cm, the 
initial solidification begins well into the hot portion of the heater 
and gradually drops down into the mid-plane of the heater as steady state 
growth is achieved. As solute is depleted near the end of the ampoule, 
the interface drops into the cold end region. In a real furnace system, 


Pace 17 



there is only a small region over which a planar isotherm, and thus 
a planar interface, will exist. A shifting interface as illustrated 
by figure 5 implies that the furnace settings must be changed in time 
to insure that the planar isotherm is also the solidification temper- 
ature if we are to obtain a planar interface. Calculations such as these 
may be coupled to thermal analysis computations to find the appropriate 
furnace settings as a function of time to achieve optimized growth 
conditions for Hg Cd Te. 


Acknowledgements — SEMTEC would like to thank Dr. S.L. Lehoczky and 
Dr. F. Szofran of MacDonald Douglas Research Laboratories for the use 
of their phase diagram and composition profile data. The work at MDRL 
is partially funded under NASA contract NAS8-33107 .. 



3.0 HIGH GRADIENT FURNACE 


3.1 HOT ZONE 

The hoc zone of the high gradient furnace will consist of a thin 
walled metal tube that is heated by passing low voltage high current 
AC power through the tube. The furnace is being designed to run in 
air or vacuum by constructing two interchangable tubes, one of tantalum 
and one of Inconel 600. The tantalum tube would require a pressure of 
10 ** torr because of the high temperature oxidation properities of 
tantalum. A tube constructed of Inconel 600 high temperature stainless 
steel will be used in vacuum and also tried in air. 

The tantalum tube assembly is a modified tantalum oversheath for 
thermowells manufactured by Thermo/Co. of Houston, Texas. The tube 
is 12'* long with flanges welded on both ends. It is produced from 
stock 5/8" I.D. .015" wall Fansteel tantalum. 

The Inconel 600 tube was constructed, from a 1" rod stock because 
of delivery lead time and size availability of the tubing. The 1" x 
12" rod stock was drilled to .75" I.D. then the wall was turned to 
0.010". A flange of 347 stainless steel 0.5 mm thick from Goodfellow 
Metals of Cambridge, England was welded to the tube ends for power 
connections and support. 

To estimate the amount of power needed to raise the tube heaters 
o 

to 1000 C only the resistivity and surface areas need be known. 


Using the following equations: 


R » PL/2 irr t 
o 

(1) 

■ 2 ir L 
t o 

(2) 

I * J efiA^ 

(3) 

R 


v - IR 

(4) 


where: * surface area 

r * outside radius 
o 

L * length 
R * resistivity 
p * wall thickness 
I = current 
e * emissivity 
5 = Stefan-Boltzman constant 

T * temperature 
v » volts 

t = wall thickness •» 

the amount of current needed for the .015" wall 5/8" I.D. tantalum 
tube heater is - 400 amps at 3.25 volts. On the other hand, the 0.010" 
wall 3/4" I.D. Inconel tube will require = 360 amps at 7.5 volts. These 
current values are in the range of available transformers and are not 


difficult to obtain. 



3.2 ADIABATIC ZONE AND GAS INLET 


The design of the adiabatic zone has been changed from that 
discussed in the second Quarterly Report. Instead of the layered metal 
strips, a ceramic block with a thin stainless steel insert in the 
center will be tested. This is designed to give a greater flexibility 
with the sizes of adiabatic zones that can be used. The ceramic will 
hopefully provide a sufficient amount of insulation between the hot 
and cold zones. 

The design of the gas inlet for the withdrawing of heat from the 
sample in the cold end has also been changed to aid in the testing of 
the high temperature furnace. The inlet has been made removable so that 
the furnace can be operated with and without the extra cooling. 

f* 

3.3 COLD END HEATER 

A modified Applied Test Systems (ATS) series 3110 laboratory 
furnace will serve as the cold end heater. The furnace has a 3/4 n I.D. 
with a 9" heated length. With the axial insulation removed from the 
ATS furnace, the tube furnace is easily mounted to the ATS furnace. 

This stackability of the ATS furnace, gas inlet, adiabatic zone, and 
the tube furnace leads to an efficient and stable arrangement. 


3.4 STATUS 


Design changes have been finalized and drawings are completed, 
{law materials have been received and machining is in its final phase. 
Assembly and system check out will be performed in the next phase of 
the high gradient furnace project. 


PROBLEMS 


There were no problems of a technical nature that would impact 
performance in this reporting period. 


PLANS FOR THE NEXT QUARTER 

In the next quarter, the results of the solution to the one 
dimensional solidification problem will be applied to derive furnace 
settings, ampoule length and sample pull rate for optimum growth 
conditions. The diffusion constant will be determined from experimental 
data. In addition, the two dimensional diffusion problem will be 
examined in an effort to couple fluid flow and diffusion into the 
interface stability problem. 

Work on the high gradient furnace will continue, with testing 
expected to be completed near the end of the next reporting period. 


SCHEDULE AND FINANCIAL PERFORMANCE 

The following pages present schedule and financial data current 
up to this quarter. Additions to and redirection of emphasis have 
modified the schedule somewhat. 


PROGRESS AS OF 20 JUNE 1980 
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Quarterly Final 

Reports Report 


FINANCIAL REPORT 


CONTRACT NAS8-33698 

The period of performance, including 
Amendment No. S/A 2 (CPFF) , is to 
November 28, 1980. 


Total Estimated Cost 

Fixed Fee 

Total Sum Allotted 


$ 57,121.00 
4,101.00 
61,222.00 


Total Direct Labor Expended 

through 6/12/80 (Voucher 19) $ 17,273.00 

Direct Materials 862.22 

Materials Overhead 86.22 

Overhead 14,336.68 

Travel 246.05 

G&A 4.778.38 

Total $ 37,582.65 


Total Estimated Cost for 
Completion through 28 Nov. 1980 

Total Amount Billed through 
6/12/80 (Voucher 19) 

Amount Received through 
6/20/80 (Vouchers 1-18) 

Fee billed through 6/20/80 


$ 19,538.35 

37,582.65 

35,147.54 

1,460.03 


Estimated progress on technical work is 65% of work completed. 
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